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ABSTRACT 

We investigate the possibility of Photospheric Radius Expansion (PRE) during mag- 
netar bursts. Identification of PRE would enable a determination of the magnetic 
Eddington limit (which depends on field strength and neutron star mass and radius), 
and shed light on the burst mechanism. To do this we model hydrostatic atmospheres 
in a strong radial magnetic field, determining both their maximum extent and pho- 
tospheric temperatures. We find that spatially-extended atmospheres cannot exist in 
such a field configuration: typical maximum extent for magnetar-strength fields is ~ 
10m (as compared to 200 km in the non-magnetic case). Achieving balance of gravita- 
tional and radiative forces over a large range of radii, which is critical to the existence of 
extended atmospheres, is rendered impossible in strong fields due to the dependence of 
opacities on temperature and field strength. We conclude that high luminosity bursts 
in magnetars do not lead to expansion and cooling of the photosphere, as in the non- 
magnetic case. We also find the maximum luminosity that can propagate through a 
hydrostatic magnetar atmosphere to be lower than previous estimates. The proximity 
and small extent of the photospheres associated with the two different polarization 
modes also calls into question the interpretation of two blackbody fits to magnetar 
burst spectra as being due to extended photospheres. 
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1 INTRODUCTION 

Photospheric Radius Expansion (PRE) events can occur 
during bursts on neutron stars when the luminosity of the 
object reaches the Eddington Luminosity, i.e. where the ra- 
diation force balances the gravitational one: 



<krGM»c 



(1) 



(where M* is the stellar mass and KTh the Thomson scat- 
tering opacity) and the large radiation pressure forces the 
atmosphere to expand outwards. For the hydrogen atmo- 
sphere of a 1.4M© neutron star, Z/Edd = 1.8 x 10 38 ergs _1 , 
while for a helium atmosphere it is twice that value. As 
a result, the photosphere moves to a much larger radius, 
corresponding to a drop in temperature T. For a neutron 
star with a modest magnetic field (up to ~ 10 12 G), Comp- 
ton scattering dominates the opacity in the atmosphere and 
various relativistic effects allow the atmosphere to expand 
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up to hundred kilometres, so that the temperature of the 
expand ed photosphere drop s out of the X-ray range alto - 
gether (|Hoffman et al.lll97g| ; IPaczvriski fc Anderson! Il986f ). 
The hallmark of PRE in neutron stars is thus a 'double- 
peaked' structure in the X-ray light curve of a burst, in which 
the flux increases to a maximum and then drops sharply 
(indicating the black-body temperature has decreased), be- 
fore rising again steeply to a slightly larger maximum as 
the bolom etric luminosity drops again and the photosphere 
contracts (|Pacz vhskil 1 1 9831 ) . 

PRE is characteristically seen in Type I X-ray bursts 
from accreting neutron stars, in which the build-up of 
accreted material leads to a thermonuclear explosion 
on the surface of the star, causing a huge increase in 
luminosity. PRE bursts have typically been used to con- 
strain the mass and radius of the neutron star, thus 
potentially co nstraining th e equa ti on of state of dense 
matter (e.g. iDamen et all Il990l: iGallowav et all |2003l: 



Ozel Giiver fc Psaltisl 12009; ISteiner. Lattimer fc Brownl 



2010l ; ISuleimanov. Poutanen fc Werne* 201ll h However, 
PRE is generically driven by high luminosities, irrespective 
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of the underlying energy source. This has led to the recent 
suggestion that it might also happen in bright bursts from 
magnetars (|Watts et alj|20foh - isolated neutron stars with 
dipole magnetic fields above ~ 10 13 G - whose bursts 
(which occur over a wide range of luminosities) are thought 
to be powered b y large-scale reconfiguration o f the decaying 
magne tic field (|Thompson fc Duncan! 1 19951 ). IWatts et al.l 
(|2010f ) argued that observing PRE in magnetar bursts could 
put interesting constraints on the emission mechanism, 
magnetic field strength and mass-radius relationship for 
magnetars. 

The suggestion that PRE might happen in magnetar 
bursts was motivated by the August 2008 observation of a 
large (L x ~ 7 x 10 39 erg s _1 ) burst from SGR 0501+4516, 
which showed a double-peaked light curve similar to those 
seen in Typ e I X- ray bursts. In their paper on this burst, 
IWatts et all |2010h laid out several criteria required for PRE 
to occur, and argued that they were in general met for mag- 
netar bursts. For PRE to occur in a neutron star, the flux 
must be emitted from an optically thick region, the radi- 
ation pressure must be sufficient to overcome gravity and 
other confining forces, the emitting region must remain op- 
tically thick during the expansion (so that the emission 
remains close to blackbody and effective temperature de- 
creases with increasing photosphere radius) , and the opacity 
must increase with distance from the star. The last point 
is slightly subtle: in order for the photosphere to expand, 
the luminosity L must remain close to the critical lumi- 
nosity (L C r) needed to balance radiation pressure with the 
confining forced]. However, both these quantities are mod- 
ified differently by the strong gravitational field, so that 
L/L CT oc 1 + z, where z is the gravitational redshift. To en- 
sure that this quantity does not decr ease with radius, which 
woul d make expansion impossible (|Paczvriski fc Anderson! 
1 1986h . the opacity (which determines L cr ) must therefore 
increase with radius. In Type I X-ray bursts, this is effected 
by the Klein-Nishina corrections which reduce the Thom- 
son cross-section at high temperatures close to the stellar 
surface. 

The presence of a magnetar-strength magnetic field 
complicates the hydrostatic expansion of the atmosphere 
in three significant ways. In the 'trapped fireb all' picture 
of magnetar bursts (|Thompson fc Duncan! Il995l ). the huge 
release of magnetic energy leads to the creation of a pair- 
plasma, so that the atmosphere is dominated by this dense 
pair-gas, rather than baryonic matter ablated from the star's 
surface. Further, closed field lines provide a strong confining 
force for both baryonic and leptonic matter, since plasma 
cannot easily move perpendicular to the field. A straightfor- 
ward calculation (see section I2.1|l demonstrates that a very 
strong field can easily confine eve n the largest giant flares 
with L ~ 10 44 erg; s _1 jLambll 19821 ). so that PRE in a mag- 
netar will likely only occur in open field line regions. 

The final effect of the magnetic field, and the most 
important one in the present work, is to modify the elec- 



tron scattering cross-sections by several orders of magni- 
tude, depending on the polarization state of the scattering 
photon. The strong magnetic field suppresses electron mo- 
tion perpendicular to the field, so that photons that try 
to excite this motion (i.e. that are polarized perpendicu- 
lar to B; the 'Extraordinary' or E-mode) have a greatly 
reduced scattering cross-section compared to the Thom- 
son scattering cross section err-h, while photons polarized 
parallel to B (the 'Ordinary', O-mode) are largely unaf- 
fected. The modified polarization-dependent cross-sections 
will increase the critical luminosity (sometimes called the 
'magn etic Eddington limit') by s everal order s of magni- 
tude (|Thompson fc Duncan! Il995l ; iMillerl 1 19951) . Addition- 
ally, since the E-mode cross-section scales roughly as (T/B) 2 
(see Equation I14jl , the opacity increases steeply with dis- 
tance from the star from the decrease in field strength, and 
becomes strongly temperature dependent. L cr is therefore a 
strong function of radius. As we will demonstrate, the strong 
temperature and field dependence of the opacity has a pro- 
found effect on the structure of the magnetar atmosphere in 
comparison to the non-magnetic case. 

The object of this paper is to explore the structure of 
a magnetar atmosp here at very high luminosities . We fol- 
low the approach of iPaczvhski fc Ander son ( 198^) and cal- 
culate the structure of a series of hydrostatic atmospheres 
with different masses, base temperatures, and magnetic field 
strengths, solving t he equations for stellar struct ure. The 
main difference from lPaczvnski fc Anderson! (| 1986T ) is in our 
consideration of the opacity, which is dominated by electron 
scattering in both cases. They consider a non-magnetic star, 
for which the electron scattering cross-section is the Thom- 
son one (modified at high temperatures by Klein-Nishina 
corrections). We instead use the cross-sections modified by 
the strong magnetic field, so that the radiation is split into 
two polarization modes, and only E-mode photons diffuse 
through the atmosphere. 

The paper runs as follows. In Sections I2.H2.2I we dis- 
cuss the concept of a 'critical luminosity' in more depth, 
and present the equations we use to calculate the structure 
of the atmosphere. Section 12.31 focuses on the electron scat- 
tering cross-sections in a super-strong magnetic field. The 
field introduces several complications (such as the depen- 
dence on photon energy of the scattering cross-section and 
the presence of the cyclotron resonance), and we explain 
how we calculate the effective opacity for a thermal photon 
distribution. Section [3] presents the main results of our cal- 
culations, and Section [4] compares our results to previous 
calculations of the critical luminosity. Finally, in Section 
we examine additional physical processes (most significantly, 
a pair plasma gas and magnetic confinement from closed 
field lines), and argue that neither of these are likely to af- 
fect the qualitative conclusions of the paper, by which we 
mean the non-existence of hydrostatic extended magnetar 
atmospheres. 



1 Here and thro u ghout the paper, we follow 
iPaczvriski fc Andersonl lll986l ) in defining L cr as the actual 
maximum luminosity as a function of radius, modified by the 
changing opacity and gravitational redshift, whereas Z/Bdd ls 
strictly given by eq[T] 



2 MODEL 

The structure of the atmosphere of a magnetar can be cal- 
culated from the equations of stellar structure. We do this 
from the surface of the magnetar to a point far away from 
the star where the atmosphere can be said to have ended. 
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We describe the various equations in detail, starting with the 
critical luminosity and magnetic field in Section ^. H and the 
stellar structure equations in Section 12.21 The form of the 
opacity in a strong magnetic field is discussed in Section [273] 
We then treat the boundary conditions and computational 
method in Section \2. 41 and briefly summarize our model. 



2.1 Critical luminosity in a super-strong magnetic 
field 

Our goal is to compute the ratio between the local luminos- 
ity L and the critical luminosity L cr throughout the atmo- 
sphere, with the latter defined as the luminosity at which 
the outward radiation force is exactly balanced by inward 
forces. If L/Lcf is greater than one at any point, the atmo- 
sphere will be unstable there, and thus not in hydrostatic 
equilibrium. 

In this work we define L CI to be the critical luminosity at 
which the radiative force exactly matches gravity, reserving 
the term Eddington luminosity for the special case where 
the opacity is given by the Thomson scattering opacity. This 
critical luminosity is given by 

-1/2 

Lcr = 11-^1 , (2) 



where i? g = 2GM/c 2 is the gravitational radius, M is the 
mass of the star, k is opacity, r is the radius measured from 
the centre of the star and c and G are the usual constants. 

In a magnetar atmosphere the opacity will be reduced 
by several orders of magnitude from the Thomson opacity 
(|Thompson fc Duncan! fl995l ; lMiHerlll995l ), giving a critical 
luminosity of order 10 40 — 10 41 erg s~ , which is in the range 
of typical luminosities of magnetar bursts, and several or- 
ders of magnitude larger than the Eddington luminosity of 
2 x 10 38 erg s _1 . This critical luminosity is inversely propor- 
tional to the opacity, which means that in a super-strong 
magnetic field it will be dependent on the magnetic field 
strength, and thus be variable throughout the atmosphere. 
We discuss the effect of the magnetic field on the opacity in 
detail in Section [2.31 

In closed field line regions, gravity will not be the domi- 
nant confining force. Here, plasma is trapped by suppressing 
the transport of charges across field lines, and so the radia- 
tion force is counteracted by magnetic stress. In a magnetar 
this effect will d ominate grav ity, giving for the critical trap- 
ping luminosity l|Lamblll982T ) 



L tI ~ 2.1 x 10 4a — 



B 
B c 



10 km 



ergs 



(3) 



where B is the magnetic field strength, R+ the radius of the 
neutron star and B cr — 4.4 x 10 13 G is the quantum criti- 
cal magnetic field strength, the field strength for which the 
cyclotron and rest mass energies of the electron are equal. 

For a typical magnetar this critical luminosity will be 
above 10 50 erg s _1 , far more than the luminosity of even the 
rare giant flares. As PRE requires the luminosity to reach 
the critical luminosity to cause expansion, it is clear that 
PRE can only ever occur in open field line regions. In the 
remainder of this work we will therefore only treat open field 
line regions. 

Treating only open field lines does not limit the appli- 
cability of our models greatly. Even in the trapped fireball 



model (jThompson fc Duncan! 1 19951 . 120011 ). where the radi- 
ation is emitted from the surface of a fireball trapped in 
the closed magnetic field lines, most of the radiation will be 
emitted through open field-line regions. This is because the 
fireball is the hottest and thus most luminous near its base, 
where the surface of the fireball borders on the open field 
line regions, so that the radiation is emitted into the open 
field line region. 

We consider a purely radial magnetic field with the field 
strength falling off as the field of a dipole directly above ei- 
ther of the magnetic poles; the modifications incurred at 
non-zero colatitudes are not expected to qualitatively al- 
ter our conclusions. We use the fully general relativistic 
equations for the radial dependence of the magnetic field 
<|Pettersonlll974l ; IWasserman fc Shapirdll983h 



B 



6/i 



R 



■In 1 



r 



+ 27 + 1 



(4) 



where p is the magnetic dipole moment. 



2.2 Stellar structure equations 

We consider a spherically symmetric atmosphere model, 
using th e general relat ivistic equations of stellar structure 
given bv lThornd (|l977h to calculate the radial structure of a 
magnetar atmosphere. The relevant equations are the equa- 
tion of hydrostatic equilibrium, the equation of energy trans- 
port, the equation for optical depth and the mass continuity 
equation. We reformulate these equations as: 



dP 
dr 



GMp 



r 



l + (4 



1-5/3)-^ 
pc z 



(5) 



as a 



specialized form of the Oppenheimer-Volkoff equation, 

(6) 
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(7) 



(8) 



where P is the pressure, T is the temperature, p is the den- 
sity, t is the optical depth, k is the opacity (Section I2.3|l 
and AM is the rest mass of the atmosphere above radius 
r. V is defined as d log Tj d log P, and depends on the 
manner in which energy is transported. In principle, this 
is a combination of radiative transport, convective trans- 
port and conductive transport. However, conduction is only 
significant compare d to radiation at much higher densitie s 
than we consider (|Potekhin. Chabrier fc Yakovlevl I2007T ). 
while convection, which is included in the models from 
IPaczvriski fc Anderso n (1986), is stro ngly suppressed by the 
magn etic field in the magnetar case (jRaiagppal fc Romanil 
Il996l ). Thus, we set V = V ra d, which is given by 



V ra d — 



16ttcGM(1 - P) 
1 + (4 -1.5/3)- 
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-1/2 
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with Loo the luminosity as seen by an observer at infin- 
ity. Assuming blackbody emission, this luminosity is linked 
to the temperature at the photosphere (where r = -R p h) 
through 



Lea = ^TvassRphTX ( 1 — 

Jtph 



(10) 



a Stefan-Boltzmann law modified by general relativity, 
which imposes an effective reduction in solid angles by the 
factor 1 — Rg/Rph- Here T p h is the temperature at r = R p h 
and <7sb the S tefan-Boltzmann constant. Th ese equations 
are those from IPaczvriski fc Anderson! (|l98rj ). with correc- 
tions to two typographical errors in Equations (3a) and (4b) 
of that work. 

We consider a purely radial magnetic field and use an 
ideal gas equation of state, so that the pressure in the at- 
mosphere is simply given by the sum of radiation pressure 
and gas pressure: 



P = P S + P T , P s = 



kTp 
M m H ' 



4ctsb t 4 
3c ' 



(11) 



where T is the temperature, p the density, p the mean molec- 
ular mass per free particle (both ions and electrons), k the 
Boltzmann constant and jtih the mass of a hydrogen atom. 
This equation for the radiation pressure assumes the gas 
is in thermodynamic equilibrium with the radiation field. 
This assumption will break down at small optical depth, 
but should be reasonable below the photosphere, which is 
the region we are most interested in. We ignore the effects 
of any mildly relativistic electrons, as these will not change 
our qualitative results. 

We assume the atmosphere to be entirely composed of 
fully ionized helium. This is not necessarily a realistic com- 
position for a magnetar atmosphere, but does make for easy 
comparison to the nonmagnetic case. We will briefly discuss 
the effect of a more realistic composition, including an elec- 
tron positron pair plasma, in Section [STT] 



2.3 Opacity 

The opacity in a super-strong magnetic field is strongly re- 
duced from the Thomson opacity, because the scattering 
cross section of one of the two polarization modes of the 
photons is suppressed by the magnetic field. Therefore, the 
opacity decreases strongly with increasing magnetic field 
strength, which means that the critical luminosity increases 
with magnetic field strength. We will now go into what this 
means for the opacity in the atmosphere of a magnetar. 

Since the photon ene rgies under consideration h ere are 
typically above 2 keV (see lWoods fc Thompsor]|2006l for typ- 
ical magnetar burst spectra), the opacity can be assumed 
to be dominated by Compt on scattering off free electrons 
|Thompson fc Duncanlll995f ). and photoelectric and atomic 
transition contributions are generally small. Some correc- 
tions for these atomic processes are needed at slightly lower 
photon energies, but these are not expected to significantly 
change our results. 

The monochromatic opacity k„ follows from the Comp- 
ton scattering cross section. In a strong magnetic field, with 
its inherently anisotropizing influence, this cross section de- 
pends strongly on the polarization mode of the photons. 
Photons in the ordinary polarization mode (O-mode or || 



mode), with their electric field vector lying in the plane 
containing their direction of propagation and the magnetic 
field, have a much larger chance of scattering than extraor- 
dinary mode (E-mode or _L mode) photons, which are po- 
larized perpendicular to this plane. The strong magnetic 
field also modifies the dielectric properties of both the atmo- 
spheric / magnetospheric plasma, and polarizes the vacuum 
|Adlerlll97lflHarding fc Laill2006l ). Accordingly, the field in- 
troduces profound modifications to the scattering opacity. 

Assuming both photon modes are polarized perpendic- 
ular to their direction of propagation, the scattering cross 
section for a photon in mode j (where j = 1 is the E-mode 
and j = 2 the O-mode) traveling at an ang le 8 with respec t 
to the magnetic field direction is given by (|Ho fc Laill2003l ) 



l + jq 



1 oj 2 (1 



Kj cos ( 



2(^- WC ) 2 + r§/4 
1uj 2 (1 + K jC os9) 2 
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+ K 2 sin 2 £ 



with 



r e = 



3m c 



(12) 



(13) 



where orh is the Thomson cross section, ui is the angu- 
lar frequency of the photon, and uq — eB/m e c is the 
electron cyclotron frequency. T c is the classical natural 
linewidth of the cyclotron resonance, and satisfies r c /cijc = 
2a f (B / B CI )/3 where otf = e 2 /he is the fine structure con- 
stant. When cue becomes a sizable fraction of m e c /ft, rel- 
ativistic corrections to this linewi dth become mandatory 
(|Baring. Gonthier fc Hardindl2005h . However, in our calcu- 
lations we use the Rosseland mean opacity, which will be 
detailed below, which means the electron cyclotron reso- 
nance gets smoothed out. The precise numerical value of 
the line width is therefore not very important, unless the 
local field substantially exceeds 10 15 G. The term Kj incor- 
porates the influences of plasma and vacuum dispersion on 
the scattering cross section, and is detailed in Appendix IA1 
where we show how this form reduces t o the various com- 
monly u s ed ap proximations, as given bv lHeroldl (|l979l ) and 
IVentural(| 19791 ). 

For radiation frequencies below the cyclotron frequency, 
this form for the cross section reduces to 



-o-CTh, 



OTh, 



(14) 



with lo the radiation frequency, cue the electron cyclotron fre- 
quency and a"Th the Thomson scattering cross section. Thus, 
the scattering cross section to O-mode photons is roughly 
constant, while the cross section to E-mode photons scales 
roughly as T 2 /B 2 when B T. The fact that the scattering 
cross sections are so different for the two different polariza- 
tion modes means that the atmosphere will have two distinct 
photospheres: one for each photon mode. As the scattering 
cross section is significantly larger for O-mode photons than 
for E-mode photons, the O-mode photosphere will be lo- 
cated outside the E-mode photosphere. 

W e use the cross sections formulated by iHo fc LaH 
|2003l ). as these incorporate the quantum electrodynamical, 
vacuum polarization and plasma dispersion effects, making 
them correct in the high-density regions of the atmosphere, 
where the commonly-used scattering cross sections in the 
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absence of dispersion from iHeroldl ll 19791 ) break down. How- 
ever, these cross sections do not incorporate Klein-Nishina 
reductions, which come into play when the electron thermal 
energy approaches its rest mass energy, i.e. kT ~ m e c 2 . This 
is also the regime where kinematic modifications to the opac- 
ity become significant, i.e., when the motions of the electrons 
must be accounted for in treating photon transport. Neither 
of these effects will be considered in the present work, as 
an opacity equation incorporating all of the above effects is 
not currently available in the literature. A side effect of this 
opacity t reatment is that it is impossi ble to reproduce the 
results of IPaczvhski fc Anderson] (jl98rj ) with our models by 



setting the magnetic field to zero, as in the nonmagnetic 
case the existence of extended atmospheres depends on the 
Klein-Nishina corrections. 

To form the total opacity, the scattering differential 
cross section in Eq. (|12|1 is averaged over angles by mul- 
tiplying it by sin 9 before integrating over 6, so as to model 
an isotropic distribution of photons. This neglects the inher- 
ent tr ansport-induced anisotropy of ph otons in the surface 
layers (jZavlin. Shibanov fc Pavlovlll995h that will introduce 
some quantitative (but not qualitative) modification to the 
outermost portions of our atmospheric profiles. 

The scattering cross sections contain two resonances: 
the electron cyclotron resonance, and the vacuum resonance. 
While these resonances are crucial for spectral modeling, we 
find that in our models they are not of great significance, be- 
ing mostly smoothed out to relatively minor features in the 
opacity due to the fact that we use the Rosseland mean opac- 
ity. As we assume the photons to be in thermal equilibrium 
with the gas throughout the atmosphere in the equation for 
radiation pressure, it is reasonable to extend this assump- 
tion and take the opacity at any point to be the Rosseland 
mean opacity at the local temperature, as given by 



1 dB v {T) 
k v dT 



Av 



dB v (T) 
dT 



Av 



(15) 



where v is the frequency of the radiation, k v is the 
monochromatic opacity and B v is the Planck function. This 
form applies when the gradients in temperature arise on 
much larger spatial scales than the gradients for opacity. As 
a consequence of averages being formed for inverse opaci- 
ties, more tenuous and low opacity domains dominate the 
determination of the Rosseland mean, so that the various 
resonance regions do not contribute much. We define the 
Rosseland mean opacity to E-mode photons as ke, and the 
Rosseland mean opacity to O-mode photons as /to- 

For the purpose of the structure of the atmosphere, the 
opacity sets the radiative energy transport through the at- 
mosphere, as described by Equation @, as well as the radia- 
tive force, which determines the critical luminosity given by 
Equation (J2J. For both energy transport and radiative force 
the relevant photon flux is the net flux, as in both cases the 
effect of identical photons propagating in opposite directions 
cancel out. Therefore, we will define an opacity K nc t, as the 
average opacity to the net outward flux of photons: 



Fq_ ( Fo 

Knet = ~^r K + I 1 — 



KE, 



(16) 



where F is the net outward flux and -Fo is the net outward 
flux in O-mode photons. This represents an opacity that is 



appropriately weighted for the relative local fluxes of the 
two polarization modes. 

In a scattering event, a photon has a certain probabil- 
ity of changing modes. These probabilities can be calculated 
from the cross sections for scattering into one specific polar- 
ization m ode and for X-rays in m agnetar atmospheres are 
given by (|Ulmerlll994l : lMilIerlll995l ) 



Fe-+o 

PO^E 



= 1/4, 



2 / 2 



(17) 



These equations are approximate and only valid in the limit 
ui <C cjc, which is true close to the star, at least up to 
a few kilometers above the surface; they also apply to the 
Thomson regime, i.e. provided uj <C m e c 2 /Tl. 

At large optical depth the net opacity will be given by 
the opacity to E-mode photons. O-mode photons will effec- 
tively be stuck in this region, as they have a much large r 
scattering cross section than E-mode photons (|Millerlll995l ). 
This does not mean there are no O-mode photons deep in 
the atmosphere, merely that they move inwards as much 
as outwards, not contributing to the net flux, and can only 
diffuse effectively by converting to the E-mode. 

At lower optical depth O-mode photons will start to 
contribute to the net flux. To quantify this contribution, we 
consider the probability an O-mode photon at optical depth 
to has of escaping the atmosphere without converting to the 
E-mode. This probability is roughly given by 



Pose = (i - Po- 



rn 



in diffusive regimes. Any O-mode photon that escapes in this 
way contributes to a net outward flux of O-mode photons at 
all points between its origin (the point where it converted 
from the E-mode) and the outer edge of the atmosphere. 

The probability given by Equation (|18[) is effectively 
zero in the region where the optical depth to E-mode pho- 
tons is greater than one, since the optical depth to O-mode 
photons is several orders of magnitude larger. Thus, as long 
as the total photon population is divided more or less equally 
over the two modes, which should be true due to detailed 
statistical balance, we can safely assume the entire net out- 
ward flux to be in E-mode photons at te = 1. Between this 
point and the outer edge of the atmosphere, these E-mode 
photons will scatter once on average, so that the fraction 
of E-mode photons that scatter in any interval drE in the 
region 1 > te > is equal to drE, assuming most of the net 
flux will still be in E-mode photons. The fraction of the net 
flux that is in O-mode photons at a point te = r' in this 
region is then given by 



Fo 
F 



Fk-s-oPesc d.TE, 



(19) 



Note that this integral is constructed from inner radii to 
outer ones, so that r' < 1 and optical depth declines with 
radius. This equation can be used to track the approximate 
fraction of the net outward flux of photons in the O-mode, 
and thus calculate K nc t. 

To test our method of determining the fraction of net 
outward photons in the O-mode we compare this to the 
analys is and Monte Carlo simulations performed by iMilleri 
(| 19951 ). The fraction of net outward photons in the O-mode 
at t = as given by Equation (|19[) should match the fraction 
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of escaping photons in the O-mode as calculated bv lMillerl . 
We have tested this for a range of values for uj/ujc and find 
that the fraction of photons that escape t he atm osphere in 
the O-mode is roughly given by 0.22a; /W [Miller! performed 
Monte Carlo simulations to calculate the radiative force di- 
vided by flux, rather than the fraction of photons escaping 
in the O-mode. However, he also performed an order of mag- 
nitude calculation to predict the results of his simulations, 
where he calculated the fraction of photons escaping in the 
O-mode, which he then multiplied with orh/c to get the 
force divided by flux. In this order of magnitude treatment 
he calculated that the fraction of photons escaping in the O- 
mode is given by O.Ioj/cjc, which is a factor two lower than 
our estimate. However, he also found that his final predic- 
tion for force divided by flux was a factor two lower than 
the results of his simulations. Thus, if we take the results 
of his Monte Carlo simulations and calculate back to the 
fraction of photons escaping in the O-mode by dividing by 
(jTh/c, wc find that our predictions almost exactly match 
his simulations. 

2.4 Boundary conditions and method 

We solve six differential equations of stellar structure, for 
pressure, temperature, atmosphere mass and optical depth 
to O-mode photons, E-mode photons and a combination of 
both, as given by Equations ([5]), ([6]) and ((Sj and three ver- 
sions of Equation Q : for the optical depth to E-mode pho- 
tons, O-mode photons, and the net flux of photons. We do 
this by integrating these equations from the stellar surface 
at Ri, = 10 km to J? cn d, which we choose as the point where 
the density becomes so small that the atmosphere has effec- 
tively ended, choosing this cutoff density as 10 -15 g cm -3 . 
This is in principle an arbitrary value, but we have verified 
that increasing this value to 10~ 12 g cm~ 3 does not change 
the results at all, indicating that our chosen value is likely 
even lower than necessary. Additionally, we find no change 
in our results from fixing R cn a at 1000 km, which is well 
beyond the expected size of the atmosphere. 

We solve Equation (|19[) for the fraction of the net flux 
contained in O-mode photons from the point te = 1 to J? en d- 
We integrate with a straightforward iterative shooting pro- 
cedure, using an eighth order Runge-Kutta meth od with em- 
bedde d Dormand-Prince error estimation from IPress et alj 
(2007). We start with guesses for all undetermined param- 
eters at the stellar surface, and iterate until all the outer 
boundary conditions are met. 

Our input parameters are the magnetic field strength at 
the surface £>,, the temperature at the surface of T* and the 
total mass of the atmosphere AMtot • Further boundary con- 
ditions are provided by the assumption that the atmosphere 
has ended at -R cn d, giving r = to — te = and AM = at 
r = -Rend- The final boundary condition is provided by the 
assumption that Fo/F = at te = 1. 

We calculate a series of atmospheres with total mass 
AMtot in the range 10 19 — 10 22 g, and surface temperatures 
from 1 keV to a few hundred keV, repeating this for mag- 
netic field strengths at the surface of the star between 10 13 
and 10 G, a range chosen to describe a typical magnetar. 
The range of masses is chosen to encompass the mass from 
iPaczvnski fc Anderson! 1 19861 ). which is ~ 2 x 10 20 g for all 
models in that work, and to span a sensible density range at 



the base of the atmosphere, ranging from 10 3 g cm~ 3 , which 
is roughly the limit for a cold neutron star, to 10 7 g cm -3 , 
which is roughly where it becomes impossible for X-rays to 
propagate through the gas, due to the plasma frequency be- 
coming greater than the typical photon frequency. 

Summarizing, our model consis ts of solving the differ- 
ential equations of stellar structure ()Thornelll977l ) from the 
surface of a magnetar out to a point where the density is 
so low that the atmosphere has effectively ended. The in- 
put parameters of our model are the surface magnetic field 
strength, surface temperature of the star and the total mass 
of the atmosphere. We make the following assumptions. 

• The magnetic field is purely radial, with the radial de- 
pendence of a dipole field right above the magnetic pole. 
This also means there is no radial magnetic pressure com- 
ponent, so that magnetic pressure can be ignored. 

• The atmosphere consists of pure fully ionized helium. 

• The radiation field is in local thermodynamic equilib- 
rium with the gas. 

• The opacity is dominated by Comp ton scatter i ng, w ith 
the scattering cross sections as given bv lHo fc Lail (|2003h . 

• Relativistic effects in scattering that can potentially 
modify the opacity somewhat for temperatures in excess of 
50 keV are neglected. 

• For purposes of the structure equations, the relevant 
opacity is the opacity that belongs to the net outward pho- 
ton flux, as given by the sum of the Rosseland mean opacities 
to E-mode and O-mode photons, weighted by the relative 
contribution of those two modes to the net outward flux. 



3 RESULTS 

Our results differ dram atically from those found for non- 
magnetic atmospheres bv lPaczvhski fc Anderson! |l98^) . We 
find no stable atmospheres with photospheric height greater 
than about 10 meter. Although the radius of the photo- 
sphere increases as surface temperature is increased, it does 
so by a very small amount. Additionally, the temperature at 
the photosphere is higher rather than lower for models with 
larger photospheric radius, because the increase in surface 
temperature required to make the atmosphere expand more 
than compensates for any temperature decrease caused by 
having a larger photospheric radius. Our results also show 
that the photospheres for E- and O-mode photons are always 
close together, both spatially and in terms of temperature, 
despite the large difference in the opacity to these two dif- 
ferent modes. An overview of input and output parameters 
for our computed hydrostatic atmospheres is given in Table 

ED 

The general structure of the atmosphere in our results 
is illustrated in Figure [T] Almost all the matter in the at- 
mosphere is contained within the first ten meters, with only 
a fraction ~ 10~ 7 of the total atmospheric mass beyond 
that region. In those first few meters, the temperature drops 
rapidly by about an order of magnitude, to stabilize from 
there on out in an unheated coronal region. The density 
rapidly drops down to the cutoff density where we define 
the atmosphere to end. The pressure is dominated by gas 
pressure in the innermost region, and becomes dominated by 
almost constant radiation pressure when the density drops. 



© 2012 RAS, MNRAS 000,[THT2] 



Models of hydrostatic magnetar atmospheres at high luminosities 7 



e o 



cfl- 3 - 



T 2 



p„/p 




-2-10 1 
log (r-R») [m] 

Figure 1. Radial structure of a magnetar atmosphere model, 
showing temperature T, density p and gas pressure P g divided 
by total pressure P plotted against radius for B = 10 14 G, 
AM = 10 20 g and fcT* =50 keV. The two vertical lines repre- 
sent the E- and O-mode photospheres. Density and temperature 
drop rapidly below the E-mode photosphere. From there on out, 
temperature stabilizes while density continues to drop, causing 
the total pressure to change from being gas pressure dominated 
to radiation pressure dominated. 



This structure is basically the same for all different combi- 
nations of T t , B and AM tot , with different values giving 
only numerical differences. In general, a higher tempera- 
ture at the surface gives a higher temperature and luminos- 
ity throughout the atmosphere, with too high temperatures 
causing the atmosphere to become unstable so that no so- 
lution to the equations can be found. Higher magnetic field 
strengths mean a higher critical luminosity due to a reduced 
opacity, so that a higher luminosity can propagate through 
the atmosphere before becoming unstable. Finally, a higher 
total atmospheric mass allows a higher surface temperature 
before becoming unstable, with a lower luminosity emerging 
from the photosphere. 

The height of the photospheric radius of a few meters 
and the density structure of our model atmospheres can be 
understood reasonably well in terms of the standard atmo- 
spheric scale height. When considered at small height com- 
pared to the stellar radius, in an atmosphere composed of 
pure helium, above a 1.4 Mq star, this scale height is given 
by 



H 



4 [keV] Cm - 



(20) 



Thus, for a typical model atmosphere with fcT* = 50 keV 
this scale height is about two meters. Over this length scale, 
the density will drop by a factor e, which causes a strong 
drop in the temperature as well, as the two gradients are 
linked in optically thick regions. However, as the tempera- 
ture drops, the scale height, which is proportional to tem- 
perature, decreases. This causes a rapidly steepening gra- 
dient in the density and temperature profiles, so that the 
density decreases by many orders of magnitude inside a few 
atmospheric scale heights, relative to that calculated at the 
surface. Over the range of this huge drop in density the ma- 
terial becomes optically thin, which causes the temperature 
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Figure 2. Stability of several model atmospheres, as given by 
the ratio between the local luminosity L (which asymptotically 
approaches Loo in Eq. (1 1 I t at large radii) and the critical lumi- 
nosity L cr plotted against the height above the stellar surface, 
for B = 10 14 G, AM = 10 20 g and fcT* = 65, 50, 30 and 10 
keV, as marked in the figure. The variations are mostly caused 
by the variations in opacity, with which L/L C r scales linearly, as 
any other radius dependent terms are just gravitational correc- 
tions. The initial downwards slope at low height is caused by the 
temperature dependence of the opacity through the photon fre- 
quency. The sudden rise between the two photospheres is where 
the net radiation force switches from being E-mode to O-mode 
dominated, causing an increase in opacity. 



profile to cease tracing the mass density profile and signifies 
the photospheric radius. 

Such gas pressure domination of the atmospheric struc- 
ture close to the surfac e of the star is also present in t he non- 
magnetic results from IPaczvriski fc Anderson! l|l986t ). How- 
ever, in their results a huge radiation pressure supported 
region, extending up to 200 km, sets up above this gas pres- 
sure supported region. In contrast, in our models, no sig- 
nificant radiation pressure supported region is present. The 
reason for this is that for a radiation pressure supported re- 
gion to be stable, the luminosity has to be almost equal to 
but below) the critical lumi nosity throughout this region 
Paczvhski fc Anderson|[l986l ) . This is also the only way to 
have an extended atmosphere, as the thermal gas would need 
an unusually high (i.e., relativistic) temperature to support 
a stable, extended atmosphere. For an extended atmosphere 
to be in hydrostatic equilibrium requires 



dP r 
dr 



= -P- 



GM 



1 _ ik 



-1/2 



(21) 



Substituting for the gradient of radiation pressure using 
Equations ((SI Inland II ip gives: 



L 



1 _ Rg 

r 



-1/2 



(22) 



This equation is of general validity, and is satisfied in the 
nonmagnetic case for a large range of altitudes for an ap- 
propriate choice of temperature due to a fine balance be- 
tween the general relativistic and Klein-Nishina corrections 
(|Paczvnski fc Andersonlll98rJ l. 

Figure [2] shows the ratio between luminosity and crit- 
ical luminosity in a magnetar atmosphere as it follows 
from our models. This figure is analogous to Figure 2 of 
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IPaczvriski fc Andersonl (|l986l ). but while their results show 
that the luminosity remains within a factor 10 -4 below the 
critical luminosity up to a height of 200 km, our results show 
a fluctuation of several orders of magnitude over that range, 
which shows that the equality given in Equation (|22l) can- 
not be met in our models. These strong fluctuations in L/L CT 
are caused by the strong variations in opacity, which are il- 
lustrated in Figure [3] This figure shows that the opacity in 
a magnetar atmosphere drops rapidly with both decreasing 
temperature and increasing magnetic field strength. 

In our models, the temperature drops rapidly up to the 
photosphere, so that the opacity also drops. The tempera- 
ture only becomes close to constant when it decouples from 
the density outside the photosphere, while for the photo- 
sphere to extend the opacity would have to be close to con- 
stant over an extended region below the photosphere. Ad- 
ditionally, at much larger length scales the decreasing mag- 
netic field strength would also prevent the opacity from be- 
ing constant over an extended range of radii. 

We thus find that extended hydrostatic magnetar atmo- 
spheres in open field line regions do not exist, as the strong 
radial dependence of the opacity prevents the existence of 
an extended radiation pressure supported region. This pro- 
vides an importa nt additional ingr edient to the results of 
our earlier work (| Watts et al.ll20ich . While the criteria set 
out in that work for PRE to occur in magnetars are all valid, 
its preliminary considerations omitted the requirement for 
near equality of the luminosity and the critical luminosity 
throughout the atmosphere. This turns out to be a criterion 
that cannot be met in open field line re gions of magnetars . 
In summation, the conditions given in IWatts et ail (|201Ch 
are indeed necessary for PRE to occur, but they are not 
sufficient. 



4 MAXIMUM LUMINOSITY 

We now consider the maximum luminosity that can prop- 
agate through the atmosphere of a magnetar as it follows 
from our results, and compare this to previous work. The 
critical luminosity in the atmosphere of a magnetar has 
been estimated several times in the literature. Two results 
in particular are gen erally quoted. The first is th e estimate 
Lcr = (ajc/cj 2 )LEdd ((Thompson fc Duncan! fl995l ). which is 
based on the reduction of the scattering cross section of E- 
mode photons by the magnetic field as compared to the non- 
magnetic case. This estimate only takes into account E-mode 
photons, based on the assumption that the entire luminos- 
ity will diffuse outwards in the form of E-mode photons, as 
these have a much lower scattering cross section. The sec- 
ond estimate we compare our results with is the one from 
iMillerl (|l995l ). which is based on the radiation force exerted 
by O-mode photons escaping the atmosphere, and is given 
by L cr = 5(wc/w)iEdd- This equation is based on a treat- 
ment in the reg ion around the O-mode photosphere, which 
is where IMillerl expects the critical luminosity to be lowest. 
Therefore, we calculate this equation at the O-mode pho- 
tosphere, and treat the result as the maximum luminosity 
to be propagated through an atmosphere without causing 
instability. Note that the critical luminosity in our models 
fully incorporates the effects of gravitational redshift, while 
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Figure 3. General behaviour of the electron and ion scattering 
opacities with differing temperature and magnetic field, illustrat- 
ing the dependencies on temperature and magnetic field strength 
of the opacity that determine the stability of our atmosphere 
models. Plotted are the Rosseland mean opacities for E-mode 
photons versus magnetic field strength, showing both electron and 
ion scattering opacities, for p = 100 g cm~ 3 and fcT* =10, 1 and 
0.1 keV, as marked in the figure. These temperatures are chosen 
to show the general temperature dependence. A different density 
would not change these curves as long as the gas is still opaque to 
X-rays, except for the position of the vacuum resonance, which is 
the small resonance peak visible in the electron scattering opacity 
around 10 12 G for 0.1 keV and around 10 13 G for 1 keV. This 
resonance is present in all curves, but smoothed out to varying 
degrees in calculating the Rosseland mean. The maximum of the 
curves lies at the cyclotron resonance, and at lower magnetic field 
than this the opacity reduces to the Thomson opacity. It is clear 
that the ion and electron scattering opacities are comparable in 
the high field limit, while the electron scattering opacity domi- 
nates at lower fields, with the crossover between these two cases 
lying at higher magnetic field strength for higher temperatures. 
In our models, the lowest temperatures are about 1 keV at 10 14 
G and about 3 keV at 10 15 G, so all our models lie in the re- 
gion where electron scattering dominates, justifying our choice to 
neglect the ion scattering opacities. 



the given estimates do not, but this does not cause major 
discrepancies. 

A comparison between these different critical luminosi- 
ties is shown in Figure 21 and shows that the critical lu- 
minosity we find is generally lower than previous predic- 
tions, although still significantly higher than non-magnetic 
limit. Our results foll ow the behavior of the relation from 
iThompson fc Duncan! up to the E-mode photosphere, but 
are an order of magnitude lower, which is caused by the fact 
that we use the Rosseland mean opacity, which takes into 
account a thermal distribution of photons, rather than the 
monochromatic opacity their estimate is based on. Between 
the E-mode and O-mode photo s phere s, our results switch to 
following the form from IMillerl (|l995T ). as here the O-mode 
photons start to dominate our opacity. 

Our results show that a single critical luminosity for the 
entire atmosphere, or the maximum luminosity that can be 
propagated through the atmosphere, is roughly given by the 
lowest luminosity out of two approximations. The first is the 
approximation in the E-mode dominated region following 
from our models, where the critical luminosity is the mini- 
mum value that 0.1(u;c /to) 2 L^dd takes on below the E-mode 
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Figure 4. Critical luminosity versus radius, comparing our nu- 
merical results to various theoretical predictions of this quantity. 
These results are for B = 10 14 G, AM = 10 20 g and fcT* =50 
keV. Note that the Miller (1995) value is a straight line, because 
Miller calculates the maximum luminosity that can be propagated 
through an atmosphere, rather than a crit ical luminosity at each 
radius. The Paczviiski & Anderson (1986) value is the nonmag- 
netic value. 

photosphere. The second is the approximation for the region 
where O-mode photons dominate the radiation force, which 
gives a critical luminosity of roughly S ucA^Eddi ca lculated 
at the O-mode photosphere, as given bv lMilIerl(|l995l ). In our 
models the former is always the lower of the two, but this 
does not necessarily have to be the case in more detailed at- 
mosphere models. These approximations break down when 
they get close to the nonmagnetic Eddington luminosity, 
where the critical luminosity will become equal to this non- 
magnetic value. 

In our results, we find that this critical luminosity for 
the whole atmosphere is never approached closely. This is be- 
cause in our models this critical luminosity is set just above 
the surface of the star, where the pressure is dominated by 
gas pressure. The critical luminosity gives the luminosity 
at which radiation force balances gravitational force, but if 
there is significant gas pressure, this will also contribute to 
the force balance. Thus, in a gas pressure dominated region 
the maximum luminosity that can be propagated without vi- 
olating hydrostatic equilibrium will actually be lower than 
the critical luminosity we calculate above, with the precise 
value depending on the atmospheric structure. 

5 ADDITIONAL PHYSICS 

Throughout this work we have made a number of simpli- 
fying assumptions. None of these assumptions qualitatively 
alters our main result: extended hydrostatic atmospheres of 
magnetars cannot exist, as near equality of the luminosity 
and the critical luminosity is not possible over an extended 
range of radii. 

Two of our main simplifying assumptions are ignor- 
ing the presence of an electron-positron pair plasma in the 
composition, and ignoring magnetic confinement effects. We 
will now briefly elaborate on why neither of these additions 
would enable extended hydrostatic atmospheres of magne- 
tars. 



5.1 Electron-positron pair plasma 

For the composition we have assumed an atmosphere com- 
posed entirely of fully ionized Helium. Changing this to a 
different atomic composition would change the results by a 
small numerical factor, and thus not have any significant 
impact on our conclusions. However, a realistic magnetar 
atmosphere will also contain an electron-positron plasma, 
particularly if it is hot enough to contain mildly-relativistic 
electrons. This could make a significant difference, as it adds 
scattering particles without adding a significant amount of 
mass, thus increasing the opacity. 

The number density of a one-dimen sional, magnetized, 
electron-positron plasma is give n by jCanuto fc Ventural 
ll977l ; lThompson fc Duncanll 19951 ) 



n c ± = 




The appearance of the crucial exponential factor encapsu- 
lates the rest mass contribution to the pair chemical poten- 
tial when the pairs are in equilibrium with photons. Cal- 
culating this number density for some of our atmosphere 
models shows that the exponential dependence on T causes 
a drop of many orders of magnitude in the pair plasma num- 
ber density in the first few meters of the atmosphere, where 
the temperature drops rapidly. This means that either the 
contribution of the pair plasma to the total number density 
of electrons in the rest of the atmosphere is negligible, or 
that the number density of electrons in the high tempera- 
ture region at the base of the atmosphere is incredibly large, 
which would then cause the critical luminosity to drop below 
the radiation luminosity. 

The rapid temperature drop in the first few meters 
of the atmosphere is also a feature of the models from 
iPaczvnski fc Anderson! (1 19861 ). and likely of any magnetar 
atmosphere model. Thus, we conclude that even including 
an electron-positron plasma in a hydrostatic model would 
not enable an extended atmosphere, although such an ex- 
tension is anticipated to influence the detailed atmospheric 
structure. 

5.2 Magnetic confinement 

The other major assumption in our model that we know 
to be wrong is treating the magnetic field as purely radial, 
which means we ignore any magnetic confinement effects. 
As we noted in Section 12.11 the confinement caused by an 
closed field line region is too strong to allow PRE, as the 
required expansion will never occur. 

In a poloidal field structure, the non-radial (i.e., con- 
fining) component of the field scales as sin O at colatitudes 
O. Even if such a small closed field component were to con- 
tribute to the confinement, Equation ([3| has the same B 2 
dependence as the critical luminosity from Compton scatter- 
ing. The same can be asserted for quasi-toroidal field compo- 
nents from non-poloidal field morphologies. Accordingly, the 
vertical gradient of the non-radial component of the mag- 
netic field would be large. Thus, just as for Compton scat- 
tering, magnetic confinement cannot create near-equality of 
the luminosity and critical luminosity. The only way for the 
field stress to keep the critical luminosity roughly constant, 
and thus equal to the luminosity, over an extended range of 
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radii, would be for the closed field component to be at most 
a very weakly-dependent function of radius. This is highly 
unlikely in any reasonable magnetic field geometry. 



6 CONCLUSIONS 

Hydrostatic extended magnetar atmospheres do not exist, 
due to the fact that the strong dependence of opacity on ra- 
dius makes the required near equality of the luminosity and 
the critical luminosity throughout the atmosphere impossi- 
ble. A hypothetical extended magnetar atmosphere will thus 
either have large regions where luminosity is greater than 
the critical luminosity, or it will have large regions where 
the luminosity is much lower than the critical luminosity, 
so that it is not supported against gravity. Therefore, such 
an atmosphere cannot exist. This is unlike the nonmagnetic 
case, where a precise balance between luminosity and criti- 
cal luminosity means the photospheric radius can extend up 
to 200 km. The fact that magnetars do not have extended 
hydrostatic atmospheres means that PRE, as envisaged pre- 
viously, cannot work, due to the fact that stable expansion 
of the atmosphere is not possible. 

Additionally, we find that the maximum luminosity that 
can be propagated through a hydrostatic magnetar atmo- 
sphere may be lower than the critical luminosity given by 
iMillerl (|l995l ). depending o n the structure of the atmosphere. 
This is due to the fact that iMiller] |l995l ) assumes the maxi- 
mum luminosity that can propagate will be set by the scat- 
tering of O-mode photons near the O-mode photosphere. 
However, in our models the maximum luminosity is set by 
E-mode photons scattering in the highest temperature re- 
gion near the surface of the star, where the scattering cross 
section is relatively large due to the high frequency of the 
thermalized photons. This means that depending on the at- 
mospheric structure, more observed magnetar bursts might 
have reached their critical luminosity than previously as- 
sumed. 

Our results have implications for interpretation of spec- 
tral fits to magnetar burst data. Magnetar bursts are typ- 
ically fit with several different spectral models, the main 
ones being two power laws, a power law plus a black body, 
two black bodies, and more recently a power law with ex- 
ponential cutoff and optically thin thermal bremsstrahlung. 
The two black body model is typically one of the best fit - 
ting of these models |Ohve et alj|2004 iFeroci et al]|2004) . 
with typical fitting parameters giving two temperatures 
around 3 and 11 keV, and typical emission region radius 
of roughly the radius of a neutron star for the colder 
black body, and an order of magnitude smaller than that 
for the hotter component, but with a large scatter in the 
sizes. While this model has mostly been presented as purely 
phenomenological, it has occasionally been suggested that 
this could be interpreted as representing the distinct sig- 
natures of the E- and O-mode photosp heres |lsrael et alj 
|200S| ; iKumar, Ibrahim fc Safi-HarrJl201Ch . This interpreta- 
tion was attractive in the sense that it opened up the pos- 
sibility of another way of measuring the critical luminosity. 
Our results show that in a hydrostatic atmosphere the tem- 
perature difference between these two photospheres is never 
more than one or two keV, and that they are very close to 
each other spatially. This makes the interpretation of the 



two blackbody model as representing the E- and O-mode 
photospheres highly unlikely. 

Our results pose the question of what does happen when 
a magnetar burst exceeds the critical luminosity. As we have 
shown, hydrostatic extended atmospheres are impossible. 
Thus, solving this problem will require dynamical, time de- 
pendent models. These models would very likely result in 
outflows, which could give rise to several forms of observable 
emission, such as the radio emission seen from the outflow 
from the SGR 1806 - 20 Giant Flare (ICameron etlH 120051 ; 
iGaensler et al.ll2005l ; iGranot et alj|2006h . or the 'Magnetar 
Wind Nebula' recen tly detected around Swift J1834.9-0846 
(|Younes et al.|[2012T ) . Any outflows would likely be very sud- 
den, due to the short sound crossing time and burst time 
scale involved. The observation that prompted this research, 
that of SGR J05 01+4516, likely rea ched or exceeded its crit- 
ical luminosity (|Watts et alj |2010| ) , and would thus be a 
prime candidate for testing any future dynamical models 
of magnetar bursts. 
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plasma dispersion and quantum electrodynamical effects, 
but not vacuum polarization. If plasma dispersion is also 
neglected, Po 1, so that K\ ~ and A'2 — > 00, i n which 
limit Equation (|12|l reduces to the form given bv iHeroldl 
(| 19791) . 

When i ncluding vacuum polarization, the parameter j3v 
is given by (|Ho fc Laill2"003T ) 
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where wc.ion is the ion cyclotron frequency, of is the fine- 
structure constant and M = ojc/^c.ion- These equations 
have been derived under the assumption m <C 1, which is 
consistent with our results, as well as v e ^ 1, which has to 
be true for radiation to be able to propagate. 



APPENDIX A: DETAILS OF THE 
SCATTERING CROSS SECTION EQUATION 

This appendix gives more detail on the term Kj from Equa- 
tion (|12|l . which incorporates the effects of vacuum polar- 
ization and plasma dispersion into the Compton scattering 
cross sections for photons in a super-strong magnetic field. 
The index j indicates the polarization mode, where j = 1 
stands for the E-mode and j = 2 for t he O-mode. This term 
is what differentiates the treatment bv lHo fc Lai (l2003h from 
older commonly used approximations (|Heroldlll979l ; IVentural 
1 19791 ). We will now detail what Kj is with and without in- 
cluding vacuum polarization and plasma dispersion effects. 
Kj is given by 



Ki=p 



l + (-l) J 1 + 
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1/2' 



(Al) 



0o 



(X>c 



0=0o0v, 
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where uup is the plasma frequency. The parameter fiy de- 
scribes the influence of vacuum polarization. If vacuum po- 
larization can be neglected, J3y — 1, and t he cross section 
reduces to the form given bv lVentural(|l979l ). which includes 
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Table Al. Properties of the computed magnetar atmospheres. Columns show the magnetic field strength, total atmosphere 
mass, surface temperature, luminosity, density at the base of the atmosphere, radius and temperature of the E-mode 
photosphere, radius and temperature of the O-mode photosphere, and range of radii where the luminosity is below the 
critical luminosity. For a given value of B and AMtot only a relatively small range of base temperatures give meaningful 
results: too low a temperature causes the density close to the surface to be too high for photons to propagate, while too high 
a temperature breaks the assumption of hydrostatic equilibrium, due to the luminosity becoming supercritical. It is therefore 
not possible to compare all different values of B and AM W at the same temperature. 
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